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SUMMARY 



The purpose of this study was die examination of the general problem 
of propellant utilization in large liquid bipropellant rocket engines, and 
consideration of the feasibility of effecting desired utilization with a particu- 
lar method of control. 

The proposed control consists of a closed loop system whereby the 
operating mixture ratio is maintained at a specified level which in turn is 
assumed equal to the loading ratio of propellants in the vehicle tanks. The 
necessary control over propellant flow to the thrust chamber is accomplished 
through the use of differentially operated control valves located in bypass lines 
about each feed pump. 

Electronic differential analyzer simulation of the application of such a 
control system to a typical large rocket engine indicates that in the low 
frequency range the overall system dynamic behavior is essentially that of 
the controller. Controller action is limited to low frequencies in order to 
avoid the possibility of interaction of flow modulation and combustion delay 
introducing unstable thrust chamber operation. 

A comparison of the use of a bypass line about each feed pump as opposed 
to a single bypass line about one feed pump indicates that the former is prefer- 
able when corrections in mixture ratio are expected to exceed 2. 5 per cent and 
It is desirable to maintain the engine performance within 1 per cent of design 
thrust A dual line control is capable of meeting this thrust performance 
requirement for corrections in mixture ratio of up to 10 per cent. 

Practical considerations of available component limitations indicate that 

ii 



a dual line bypass control system is capable of response times in the order 
of . 07 to . 11 seconds and an expected accuracy of . 5 per cent. 
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SYMBOLS 



area 

thrust, pounds 

specific impulse, seconds 

constant 

characteristic length, feet 
pump speed 

pressure, pound per square inch 
flow rate, cubic feet per second 

loading ratio, pounds of oxidizer per pound of fuel in rocket tanks. 

universal gas constant 

absolute temperature, degrees Rankine 

volume, cubic feet 

weight, pounds 

transfer function 

molecular weight, pounds per mole 

actuator time constant, seconds 

characteristic velocity, feet per second 

gravitational constant, 32. 2 feet per second per second. 

altitude, feet 

constant 

length, feet 

mass, slugs 



p 



r 

s 

h, 

*c 

*8 

A 

t 

p 
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Cl) 

7 



d 

at 

mixture ratio 

d_ 

dZ 

time of burn out, seconds 
combustion delay time, seconds 
chamber stay time, seconds 
increment of change 
ratio of specific heats 
error signal 

density, slugs per cubic foot 
valve position 

angular frequency, radians per second 
damping factor 



Subscripts 
b bypass 

c chamber 

f fuel 

i injector or initial value 

o oxidizer 
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A PROPOSAL FOR THE CONTROL OF 



PROPELLANT UTILIZATION IN A LIQUID 
BIPROPELLANT ROCKET ENGINE 

INTRODUCTION 

The performance of a rocket propelled vehicle is greatly influenced by 
its mass at any given instant of time. The advent of large rockets, where 
structural weight is but a fraction of the propellant weight, coupled with 
strict performance tolerances requires that the utilization of the propellants 
be subject to extremely close control. 

Early attempts at establishing this control consisted primarily of cali- 
bration of the engine by ground testing. This method, while able to compensate 
for various component tolerances. Is unable to correct for any in-flight dis- 
turbances such as acceleration effects or large changes in ambient conditions. 

The purpose of this study was to investigate the use of a closed loop 
control system in order to maintain the propellant utilization at some desired 
level. The proposed control system is based on the use of bypass flow control 
around each propellant feed pump. The application of the proposed control 
system to a typical rocket propelled vehicle was studied using, in part, analog 
computer simulation. 

Rocket engine simulation has previously been treated at length by Gore 
and Carroll, (Ref. 1.); B. N. Smith, (Ref. 2.); and Sanders, Novich, and Hart, 
(Ref. 3. ). The methods used in this study represent a modification of those 
presented by Smith. The theory used to describe rocket engine performance 
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as a function of propellant flow and mixture ratio represents a simplification 
of similar considerations by Crocco, (Ref. 4. ). 

This study was undertaken by T. M. Kastner and C. D. McCullough as 
a joint thesis project to satisfy, in part, the requirements for a Master of 
Science Degree in Aeronautical Engineering at the University of Michigan, 
Ann Arbor, Michigan, September 1956 - June 1957. 

THE NATURE OF THE PROBLEM 



OBJECTIVES 

The nature of any propellant utilization system depends on the objective 

of its inclusion as part of die rocket engine. This objective, while not 

singular in concept, is determined in the main by the performance requirements 

of the rocket vehicle. The overall objective may be die attainment of one or 

< 

any combination of the following requirements: 

1. Maintaining engine operation within certain specified limits. 

2. Continuous control of vehicle mass. 

3. Obtaining the maximum energy transference from the propellant on 
board. This will call for total propellant utilization during the 
flight. In the case of a bipropellant rocket this may be interpreted 
to mean that both fuel and oxidizer are exhausted simultaneously 
with engine bum out. 

4. The future propellant requirements on a given flight may establish 
a need for varying degrees of propellant utilization. 

It might seem that these requirements can equally well be met by the 
proper choice of design parameters in the basic engine system. Unfortunately 
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the absolute achievement of any performance level is an impossibility. In 
reality there are only degrees of achievement within certain error limits. In 
order to better understand the problem of propellant utilization it is necessary 
to examine the nature of these errors and their effect upon engine and vehicle 
performance. 

CAUSES FOR DEVIATIONS FROM PERFECT PROPELLANT UTILIZATION 

In order to effect the control of propellant utilization it is necessary to 

consider the causes for the pertubations from the original engine design 

conditions. Improper propellant utilization is most evident in a shift from 

the design mixture ratio where operation at other than the design mixture 

ratio may be attributed to the following: 

1. Variation in propellant density. The basic engine design assumes 
specific values for the densities of the oxidizer and fuel. These 
densities are based upon assumed flow temperatures. If the pro- 
pellant temperatures differ from the assumed values there may be 
a change in the operating mixture ratio. This can be shown by 
considering the flow at the thrust chamber injector where 




and 




The mixture ratio may be expressed as 




or in differential form 
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Considering some variation in propellant temperatures AZ and 
AT f . with 




the change in r resulting from the temperature induced density shift 
may be expressed as 
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If an identical temperature shift is assumed for each propellant 




or in terms of some perturbation from the design value T, 
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From the above it is easily seen that for any condition other than 
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a shift in mixture ratio will result from a change in propellant 
injection temperatures. This effect is shown in Figure 1 for several 
different propellant combinations. Note that only in the case of nitric 
acid and aniline is this condition approached. It may easily be shown 
that for a LGX-ETHYLAC combination 



~ = 353 « /O ' 5 /<, 

— - .35^ 10 /o R 
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Here a 50° temperature shift will result in a 1. 3% change in r. It 
is surmised that some statistical variation in component densities 
may exist depending mainly upon supply conditions. 

2. Variations in initial loading ratio R. In general the initial loading 

of the rocket will vary statistically with the probability of any specific 
loading ratio given as 



where R is some increment of measure and g(R) is the distribution 
function for R. The nature of g(R) will depend on the accuracy with 
which the initial propellant loading can be determined. This particular 
perturbation would seem to lie outside the realm of total correction 
by propellant utilization system. However, if the nature of g(R) is 
known it is possible to reduce the likelihood of any resulting error to 
a minimum by controlling the operating mixture ratio such that there 
is an equal probability of having an excess of either fuel or oxidizer 
at bum out. This consideration is dealt with at greater length by 
Brousseau in Ref. 5. 

3. A similar statistical perturbation will result from production toler- 
ances in the system components. Here the probability of any operating 
mixture ratio will be of the form 



where g(r) is the distribution function of r. It is possible to correct 
this perturbation within the limits of accuracy of measurement of the 
instanteous of operating mixture ratio and the capability of the control 
system. Indeed, one principal advantage of any such control system 
is the fact that its use might allow relaxation of production tolerances 
of certain engine system components and thus enhance over-all production. 

4. Accelerations experienced by the rocket during its flight One effect 
of acceleration is such as to vary the hydrostatic contribution to 
injection pressure of each propellant component Here the resulting 
change in mixture ratio can be shown to be approximately 












Most of this effect is due to variations in the inlet pressure of the 
feed pumps due to the change of hydrostatic head in the propellant 
tanks. The effect may be compensated by adjusting the tankage and 
line size in such a way to insure that 
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Reichel in Ref. 6. considers the advantages of a concentric tank 
arrangement towards this end. Other possible schemes include a 
variable tank pressurization system to insure constant pump suction 
head or use of special pump design and compensation. 



EFFECTS OF IMPROPER UTILIZATION 

Justification for the use of a propellant utilization control system may be 
best effected by considering the results of improper utilization. If a rocket 
vehicle, such as that of Appendix I, is initially loaded with propellants ap- 
portioned so that engine operation at the optimum mixture ratio will result 
in complete burning of all propellants during the flight any improper utilization, 
as represented by operation at other then the designed mixture ratio, will 
result in the following: 

1. There will be residual propellant on board at burn out which will 
act to decrease the effective mass ratio of the rocket and thus its 
performance. 

2. The burning time will be decreased. 

3. There will be a deviation from the design specific impulse which 
will change the performance. 

4. The shift in mixture ratio, if large, may cause temperature limits 
of the engine to be exceeded. 

5. Improper propellant utilization may result in low frequency com- 
bustion instability as evidenced in the experimental work Barres 
and Moutet,(Ref. 7). 

6. Any variation in the total mass flow rate of the propellant will affect 

the engine thrust. This effect is generally of minor importance and 
can either increase or decrease performance. 
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Appendix n demonstrates some of the more marked effects of improper 
propellant utilization on a typical large rocket. Other examples are given by 
Rosen,(Ref. 8. ) and Reichel, (Ref. 6. ). In Ref. 8, Rosen vividly Illustrates the 
performance problem by calculating the expected performance of three theoretical 
satellite proposals with an allowance of . 025 of the initial propellant weight 
remaining on board at burn out. (. 025 is an average figure for the initial 
Viking series. ) The corrected orbits fell far short of the ideal calculations 
based upon perfect utilization and in one instance orbital velocity was not 
attained. 



PROPOSED CONTROL SYSTEM 



DESCRIPTION 

In conducting this study it was necessary to assume requirements that 

seemed typical to present applications of the large rocket engine. The essential 

objective of the control system proposed is the ability to effectively control the 

propellant consumption throughout the engine burning time in such a manner as 

to insure complete utilization of all propellant components at burn out. Other 

considerations are that engine performance with respect to thrust and chamber 

« 

pressure should be essentially constant and that the system be compatible with 
stable engine operation throughout its control range. 

Any such system must of course be restricted in its application to instal- 
lations in which the added weight and complexity of the system do not offset the 
advantages gained by its employment. 

There are several means by which control of propellant utilization may be 
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effected. Basically these consist of first, some means of measuring any 
departure from the desired utilization criteria and second, some means of 
compensation for the detected error. The measurement problem may be 
considered as one of determining the propellant on board and of determining 
the instanteous operating mixture ratio. Use of flow rate measurement 
methods seems to offer the most hope of accomplishment with the least error. 
Compensation for improper utilization will generally require that the propellant 
in excess be disposed of in order to avoid the performance limitations imposed 
by its addition to die rocket weight This may be accomplished by jettisoning 
or by burning in the thrust chamber. Jettisoning of fuel from a rocket with the 
engine operating is generally not desirable. It means a waste of fuel energy 
and, it does not correct any deviations in operating mixture ratio. It does 
have the advantage of not disturbing the thrust chamber flow. The second 
method is basically one of maintaining the mixture ratio at a specified level 
by modulation of the flow rates. A reasonably rapid response in such a system 
will avoid any variation in desired utilization provided the specified operating 
mixture ratio and the loading ratio of the rocket are identical. This assumes 
a knowledge of the loading ratio at time of lift off through launch pad instru- 
mentation and a corresponding setting of the operating mixture ratio. 

The engine and control system as shown schematically in Figure 2. In 
this system the feed pump speeds are assumed to remain constant. The flow 
rate in each propellant line at the injector is varied by using a controlled 
bleed or bypass from the injector inlet to the pump suction. Mass flows are 
obtained from the volumetric flow rate output of vane- type flow meters which 
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are calibrated against propellant temperature, 



where 



m = f><? 

/o= -P(t) 



The method of controlling bypass flow is shown in Figure 3. Here the 
operating mixture ratio is compared with the desired ratio, R, and the differ- 
ence or error signal is used to position the flow valves. The valve motions 

of the control bypass lines are opposite in sense. This differential action is 
« 

made a function of the error signal by inversion of the signal supplied to one 
valve actuator. Another possible solution is the use of a single actuator with 
differential gearing of the valve action. 

In determining the nature of the control system components several 
assumptions regarding the effects of the control system upon the response of 
the engine are made. 

1 . Controlled flow rate perturbations are to be low enough in frequency 
so that any phase lag occuring due to combustion delay time will be 
negligible. The delay time is assumed to be approximately 0. 002 
seconds, (Ref. 9. ). Therefore, it is required there be no control 
perturbations with a frequency much greater than 50 radians per 
second. 

2. It is desired that the thrust chamber pressure not vary more than 2% 
from the design value. 

3. The control system should be capable of correcting for possible 
deviations of 10% in operating mixture ratio. 

4. The rate of change of mixture ratio correction is to be limited such 
that there is no more than 1% variation in mixture ratio between injector 
flow and chamber exit flow as a result of control action. This requires 
an error signal limitation of 10% "r„ Also, assuming that the maximum 
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corrective action is demanded in the form of a step response, if the 
rate of change of mixture ratio is held so that there is no more than 
a 10% response to the error signal within a period equal to the 
chamber stay time the rate limitation will not be exceeded. This 
allows the assumption of uniform conditions throughout the chamber. 

The limitations regarding frequency and rate of flow modification as a 
result of controller action are necessitated by stability considerations of the 
basic system composed of thrust chamber and feed lines. Studies of the low 
frequency dynamics of rocket engines show the dependence of stability upon 
the phase lag between injector flow and chamber pressure. (See References 
7, 9, 10, 11, and 12). The contribution of combustion delay time to this phase 
lag is of primary importance in stability considerations. In as much as the 
behavior of the delay time under conditions of varying mixture ratio as well 
as overall mass flow is not yet fully understood it was considered desirable 
to restrict the control flow perturbations such that the phase lag due to any 
delay time would be less than one tenth radian. 

SYSTEM EQUATIONS 

Derivation of the equations describing the system is based on linearization 
about the design operating conditions. The following assumptions are made. 

1. System components may be described by lumped parameters. 

2. Suction heads for both pumps are constant 

3. System components are inelastic. 

The equations describing the system are shown in Figure 4. 

Pumps : The feed pumps are centrifugal flow pumps with typical head 
versus flow characteristics. Stepanoff, (Ref . 13) gives the following analytical 



expression for this type pump characteristic: 
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p 0 -p s = Ai'J* +8QM +CQ^ 

Given constant N and P e the pump behavior may be expressed 

AP 0 = f ! AQ K P AO 

Similar characteristics were assumed for both fuel and oxidizer pumps. 

Main Propellant Line : For small perturbations in flow the effect of line 
friction may be neglected with the only pressure drop due to the flow inertia. 
Thus 

Po-n - (&♦<?.) 

or 

APo-AP. = K, (AQ. + AQ k ) 

Injector : Injector flow is given by the equation for orfice flow as 

9 . 

or in linearized form as 

AP x -AP t = - 0 — AQ, = < k AQ. 

Flow Sensing: The problem of accurate flow sensing, or measurement, is 
one of great importance in effecting any satisfactory propellant utilization 
control system. References 13, 14, and 15, describe several methods of 
approach. The dynamics of a vane type sensing instrument are considered in 
detail by Grey in Reference 14. Time constants of, from . 5 to 3 milliseconds 
have been determined for such instruments. Therefore, the dynamics of the 
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flow sensors are assumed negligible with respect to the overall system 
dynamics. 

Bypass Line: With similar assumptions as those for the main line 

AF3' * 

Bypass Valve : In terms of valve position and flow the pressure drop across 
the valves is expressed as 



where 0 represents valve position. In perturbation form 



- 2 lfLB.'^A6 = AOb - Ko A© 



Oi. e 

Thrust Chamber: The linearized equation representing the combustion 
chamber is developed at length in Appendix in and is given as 



(KS)A^ -- K^AQ^-HC* AQ.t +K c .(l+2s)Ar 



Mixture Ratio : The mixture ratio expressed in linearized form is 

Ar * ~ AQu - AO.t 

Application of these equations to both propellant lines and the elimination 
of the pressure terms except for P c results in the following expression for the 



system: 
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The constants K L etc. refer to a change of time variable T- -p to corre- 
spond to the time scaling of the thrust chamber equation. 




While it is possible to use wholly analytical means to solve for the dynamics 
of the system, the solution of the determinants involved is laborious. Therefore, 
the use of analog computer simulation was chosen. It was considered desirable 
to check the solution for operation at design mixture ratio using the matrix 
equations to determine an equivalent transfer function representation of the 
mixture ratio response to a change in control valve position. This transfer 
function is of the general form 






A Bode plot of the frequency response for the ideal case of r r i.V] is shown in 
Figure 5. It is of interest to note that over the range of frequencies for which 
the control system is desired to act the transfer function is effectively constant 
and the phase shift is small. 



* For dual valve control , 



A© = AS. ■=. —A©* 
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The control valve proposed Is a gate type valve with an electric servo 
drive. A schematic of the valve is shown in Figure 6. A triangular orifice 
was chosen to insure constant discharge coefficient over the range of valve 
operation. Denoting the valve position for design flow condition, r = 1.37, 
as 7, valve positioning by controller action is represented in terms of non- 
dimensional valve movement . In order to obtain linearity between 

valve orifice area and actuator position, it is assumed that a nonlinear link- 
age exists such that 

A A _ A& 

A Q 



The over all valve and actuator motion may then be expressed as 

as 

6 p(Tp+0 

HereT represents the ratio of valve and actuator inertia to damping. K^. is 
the controller gain, and e r is the error in operating mixture ratio. In pertur- 
bation form 

6, = ar- 



Scaled to the same time, 




as the previous equations. 



A0 « 

e 



Kvt 

S(T r s+l) 



where 



Kv* * KV t <, 




When dual valve control with differential valve action is used 
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SIMULATION 

The computer simulation of the control system operation consisted of 
several phases, namely: 

1 . Study of die open loop dynamics of the basic system less controller 
over the expected range of operation, 

yi U 

2. Observation of the variation in equilibrium conditions as a result 
of control valve operation. 

3. Observation of the closed loop response of the system for a step 
input of desired mixture ratio. 

The computer circuit used to simulate the basic system is shown 
schematically in Figure 7. The basis for the computer circuit is the previ- 
ously developed set of linear perturbation equations. The value of all con- 
stants except K^, K^, and K c * remain essentially constant over the range of 
flow conditions considered and were thus allowed to remain as calculated for 
operation at ideal design conditions. The values of Kg, K^, and Kg* were 
varied according to the equilibrium conditions about which the perturbations 
were assumed. 

« 

In the first phase a sinusoidal perturbation in — =- was applied to obtain 

9 

the frequency response of the feed lines and thrust chamber. This was done 
about equilibrium conditions for r - 1. 37, 1. 50 and 1. 26 which represented the 
optimum and extremes of operating limits. 

The second phase consisted of the determination of equilibrium conditions 



existing over the proposed range of control action. This was done by a 
summation of perturbation effects. Starting at the design condition, successive 
increments of A r . . «t were used as a basis of calculation. This procedure 
was done for the following three modes of operation: 

1. Bypass control in each propellant line 

2. Bypass control in the fuel line only 

3. Bypass control in the oxidizer line only 

The final phase consisted of selecting a controller function which would 
give optimum response within the flow rate limitations dictated by the initial 
assumptions. This was done in terms of the system response to a step input 
for assumed values of valve actuator time constant T. of . 01 and . 02 sec. 
in addition, the frequency response of the closed loop system was obtained 
in order to observe the attenuation in flow variation with error signal frequen- 
cy. The computer representation of the valve actuator and controller Is 

shown in Figure 8. 



RESULTS AND DISCUSSION 

The results of the computer simulation for the first phase are shown in 
Figure 9. There was no attenuation in response for frequencies less than 
180 radians per second. At this point the attenuation characteristics exhibited 
a -6 db. per octave roll off with no appreciable peaking at the break point. 

This attenuation characteristic was practically identical for both the extreme 
and optimum mixture ratio conditions. The accompanying phase lag of the 
mixture ratio response is small and approximately linear over the desired 
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frequency range, oji 50. Also shown in Figure 9 is the response predicted 
by analytical calculations for r - 1.37. Lack of complete correlation may 
be ascribed to measurement inaccuracies and simulation simplifications. 

These results lead to the initial conclusion that the basic engine design is 
stable over the desired range of operation, as predicted by analytical consider- 
ations. The desired control operation frequencies lie well below the break 
frequency of the line-chamber response and the phase lag is almost negligible 
in this region. From this it is evident that within the frequency range of 0 
to 50 radians per second the response of the control system will be governed 
essentially by the response of the controller and valve actuator. 

The results of the second phase of the computer simulation are shown in 
Figures 10, 11, and 12. Of primary interest is a comparison of the results 
for dual valve control with those for single valve control. As might be expected 
dual valve control offers a wider corrective range in mixture ratio for a given 
range of valve positions. In addition it should be noted that the dual valve 
control gain as represented by the slope of the mixture ratio versus valve 
position curve is not only larger but also more nearly constant for all valve 
positions than is true for single valve operation. In single valve operation there 
is a tendency for the gain to drop as the valve is opened. With dual valve con- 
trol this effect is balance by the differential valve action. Observing the chamber 
pressure for each type of control it is seen that for dual valve control there is 
a gradual pressure drop for corrective action in either direction which is at most 
only five psi. at the extreme limits of the valve position. The control system 
requirements call for a chamber pressure variation limit of 2%. The variation 



18 



for dual valve control is well within this limit. Original system concepts 
had included a chamber pressure feed back in which the dual valves would be 
actuated in unison as well as differentially to compensate for any large vari- 
ations in chamber pressure. In view of the small variations occuring within 
any feedback provisions it is felt that over the correction range chosen this 
refinement is not necessary. The differential action of the valves keeps the 
mass flow into the chamber nearly constant and most of the pressure drop 
experienced is the effect of shifting the mixture ratio. On the other hand, the 
single valve control produces a more marked effect on chamber pressure, due 
primarily to the change in total mass flow necessary to produce a change in 
mixture ratio when only the oxidizer of the fuel flow is controlled. Due to 
greater propellant density this effect is more pronounced for a single valve 
control in the oxidizer line. The curves of flow rates in the various lines of 
the system clearly show the large bypass flow rates associated with single 
valve control where the reduced gain results in the necessity for greater 
valve movement for corrective action. It should be noted that the bypass 
flows vary linearly with valve position as might be expected from the minor 
variation in the pressure drop across the valves. Also of some interest 
is the variation in the uncontrolled flow where a single control valve is used. 
Here the chamber pressure shift caused primarily by the variation of flow 
in the controlled line results in a slight change in injection overpressure for 
the uncontrolled line and this a small change in flow. This variation is always 
in such a direction as to aid the desired mixture ratio shift but the effect is 
quite small. Using the analog simulator data for the dual valve control 
system, the variation in engine thrust, F, from the design value of 150, 000 
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pounds at r = 1. 37, was calculated for the operating range of mixture ratios. 



F = hp™* 

In the range, 1. 25 < r < 1.47, the variation was negligible. The maximum 
variation was found to be 0. 7% at r = 1.5. 

As previously mentioned the object of the final phase was to select con- 
troller parameters that would produce optimum system response commensurate 
with the desired requirements. The requirement limits are: 

1. The control system be relatively insensitive to input frequencies 
greater than 50 radians/ second. 

2. The maximum response rate to a step error signal of amplitude A, 
should not exceed A/10 Over a period of . 004 seconds. 

Examination of the system diagram shown in Figure 8 indicates that if 
the chamber and feed lines are considered as a constant and only proportional 
amplification of the error signal is used the overall system is second order 
where 



- 



s(T t $*'Vk* 






with a natural frequency of 




td* = 



UJn* 

TT 



and damping ratio of 

J 

Here the time constant T^is that of the valve actuator while the gain Ky 
is composed of both actuator gain and basic system gain. 

In order to meet the first requirement it was desired to select a valve 
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actuator with a natural frequency of approximately 50 radians per second. 

This insures an output attenuation of -12 decibels per octave beyond this point. 

To meet the second requirement it was considered desirable to limit the 
response rate such that the maximum rate would be no more than 10% of the 
amplitude of a step input over the stay time of the chamber, . 004 seconds. 

The final criterion applied to desired response was that of best possible 
settling time; that is, the time for any error to be reduced to within 5% of 
its initial value. 

Computer solutions for system response to a step input are shown in 
Figure 13 and 14, for several values of K y and T. Figure 15 gives the vari- 
ations of the response parameters as a function of controller gain for the two 
assumed actuator time constants of . 01 and . 02 seconds. The most notice- 
able variation between the two conditions is the decrease in settling time 
associated with the smaller time constant Use of a valve actuator time 
constant of . 01 seconds allows the realization of a setting time of . 068 
seconds as compared with the optimum time of . 108 seconds for the case 
of an actuator with T = . 02 seconds. In each case the controller gain for 
the optimum response is approximately the same, Ky = 250. The natural 
ability of the control mechanism to filter out inputs with frequencies greater 
than 50 radians per second is less for the faster responding control. Figure 16 
shows the closed loop frequency response of the system with control parameters 
adjusted for optimum response with an actuator time constant T = .01 seconds. 
The frequency response attenuation falls off at -12 decibels per octave from a 
break frequency of 70 radians per second rather than the desired value of 50 
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radians per second. This is correctable by using a lower controller gain at 
the expense of a greater settling time; use of a valve actuator with a larger 
time constant, again at the expense of settling time; or use of a network filter 
to operate on the error signal. This last method presupposed complete actuator 
response up to 50 radians per second. For both values of the actuator time 
constant rate limitations were met when the settling time was a minimum . 

Observation of the system response for initial conditions other than those 
at design mixture ratio of 1. 37 revealed little change in the response charac- 
teristics. The only noticeable change was a small increase in damping due 
to the slight drop in the basic system gain associated with increasing valve 
actuator position. 

The possible use of a dual bypass control system to effect stabilization 
of low frequency chamber pressure oscillations as well as performing as a 
propellant utilization control may be inferred from the study by Marble and 
Cox, Reference 10, where the use of a variable feed line capacitance was 
proposed. Bypass control is analogous to such a system if one considers that 
• • • c/C 

To be effective such a control would require good valve actuation response 
up to frequencies of 80 cycles per second, which calls for actuator time 
constants in the region of T = . 0015 seconds. This is one order of magnitude 
beyond the assumptions made in this study. 

Of considerable practical interest are certain "fail-safe" features of the 
proposed control system. In the event .of a valve or actuator malfunction the 
maximum resulting change in propellant flow will depend upon the area of the 
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control valve orifive. For the application of such a system to the rocket vehicle 
of Appendix I the orifice areas required for operation at ideal conditions were: 

A^ = . 266 in. 2 , X = . 303 in. 2 

These values correspond to valve positions of 6 and vary directly as 9. 
Examination of Figures 9, 10, and 11 indicates the feasibility of limiting 
the result of complete valve failure in the full open position by the choice 
of the available orifice area. Limitations on results of a complete closure 
may be obtained through proper choice of pump characteristics. Should 
only a single valve malfunction the control system will still be capable of 
limited corrective action through the use of the remaining valve. Provided 
that there has not been a complete failure, the remaining valve may well 
be able to both compensate for the effects of the failed valve and also 
effect a limited degree of propellant utilization. 

CONCLUSIONS 

From the preliminary investigation it is apparent that a propellant 
utilization scheme based on varying the operating mixture ratio is feasible. 

The allowable range over which the mixture ratio may be varied is a function 
of the type of propellant combination used in the engine. 

Based on linearized simulation study of the proposed system as applied 
to a typical large rocket using a LOX-ETHYLAC propellant combination, 
the following conclusions are made: 

1. Where the controller response was specified to avoid appreciable 
phase lags between line and chamber flows due to the combustion 
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time delay, the dynamics of the combined engine and propellant 
utilization system were determined, in the main, by the response 
characteristics of the controller. 

2. Within the accuracy which the loading ratio is known the system is 
capable of producing empty burn out 

3. The "best response" design requires a controller gain of Ky - 250, 
with a time constant of T e .01 seconds. This response exhibited 
no overshoot and a settling time of . 068 seconds. 

4. Dual valve control is preferable to single valve control due to the 
fact it allows a larger range over which the mixture ratio can be 
varied without exceeding chamber pressure limitations or appreciably 
effecting thrust performance. 

5. Engine performance in terms of thrust remains essentially unchanged. 
The maximum variation in thrust performance over the proposed range 
of mixture ratio control was 0. 7%. In general the actual variation was 
negligible except at the extremes of the allowable range. 

6. The system operation does not result in undesirable variations in 
chamber pressure. The maximum variation in chamber pressure over 
the range o>f mixture ratio control was 1% of the design value. 

7. With a chamber pressure feedback in the control loop for dual valve 
control the allowable range of mixture ratio excursions may be 
extended. 

8. If controller response times in the order of . 002 seconds can be 
made practicable it appears the bypass control system for propellant 
utilization could serve the additional function of controlling low frequency 
combustion instability. An immediate consequence of this would be a 
desirable reduction in injection over-pressure. 

9. The proposed system of dual valve control has a reliability advantage 
over single valve control is that corrective action may still be ac- 
complished despite the malfunction of a single valve. 
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PROPELLANT COMBINATION 




Lox and Avgas 
Lox and E thy lac 
H2O2 and Avgas 
RFNA and Aniline 
WFNA and Aniline 



1.5 

1.3 

1,82 

1.56 

1.49 



FIGURE 1 



Temperature gradient ratios for several 
typical liquid propellant combinations. 
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FIGURE 2 Schematic diagram of bipropellant rocket engine with bypass flow controls. 
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FIGURE 3 Schematic diagram for bypass flow control system. 
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FIGURE 4 System equations for propellant utilization control. 
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FIGURE 5 Predicted frequency response of mixture ratio to bypass valve 
position, open loop 



6 y V/UWf Post rtorT 



= n 



e--e 



e--o 



t 




FIGURE 6 



Schematic diagram of bypass valve 



PHASE 



30 




* 

* 

o 

* 

<0 

K 

$ 



t 

8 



K > 



s 

I 



FIGURE 7 Analog simulator diagram for propellant utilization system without controller dynamics 
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FIGURE 8 Block diagram and simulator diagram for bypass valve control 
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FIGURE 9 



Actual frequency response of mixture ratio to bypass 
valve position, open loop 
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FIGURE 10 



Steady state system operation for dual valve control 
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FIGURE 11 Steady state system operation for bypass control in oxidizer 
line only 
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FIGURE 12 Steady state system operation for bypass control in fuel line only 
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FIGURE 13 Step response of mixture ratio for different values of 
gain; ^ z . 01 seconds 



37 




FIGURE 14 Step response of mixture ratio for different values of 
gain; = . 02 seconds 
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FIGURE 15 Effect of controller gain on response parameters 
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FIGURE 16 Frequency response of mixture ratio for entire system 
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APPENDIX I. 



NOMINAL DESIGN VALUES FOR A TYPICAL LARGE ROCKET VEHICLE 

Thrust, lb. 150, 000 

Take-off thrust- weight ratio 1. 5 

Empty weight, lb. 20, 000 

Take-off weight, lb. 100, 000 

Propellants, Liquid oxygen-ethyl alcohol (with 25% HgO) 

Flow rates at optimum mixture ratio, lb. /sec. 

Main oxidizer line 375 

Oxidizer bypass line 28 

Main fuel line 274 

Fuel bypass line 21 

2 

Cross-section areas, ft. 

Main oxidizer line . 2430 

Oxidizer bypass line • 0195 

Main fuel line . 2340 

« 

Fuel bypass line . 0190 

Line lengths, ft. 

Pump to injector 4 

Bypass line 6 

Fluid velocities in propellant lines, ft. /sec. 20 

Alcohol density, slugs/ ft. ** 1.675 

Q 

Liquid oxygen density, slugs/ft. 2.210 
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Injection overpressure, psl 100 

Thrust chamber parameters 

Characteristic velocity, ft. /sec. 5,600 

Characteristic length, ft. 8 

Chamber pressure, psi 500 

Optimum mixture ratio 1.37 

Specific impulse, sec. 250 

Chamber stay time, sec. .004 

2 

Throat area, in. 209 

3 

Chamber volume, ft. 11. 6 
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EFFECT OF PROPELLANT UTILIZATION ON PERFORMANCE 

In many cases, depending on the application, it is desirable for a rocket 
to arrive at the burn out point in a minimum weight (mass) condition. For the 
bipropellant liquid rocket this optimum condition is achieved if both propellant 
components are exhausted at the same instant. Any deviation from this optimum 
will manifest itself in the form of increased burn out weight due to the mass of 
the excess component on board where the penalty of the increased burn out 
mass is a loss in performance. 

Ideally, if perfect control were maintained over the launch loading ratio, 

R, (where R is the ratio of the initial oxidizer weight to the initial fuel weight) 
and the burning mixture ratio, r, an empty burn out should be achieved. If 
the initial weight of the propellant is W p then 



w f -W>i + Wf. - VvC (H id- W fi (/4R) 



and the burning rate of the propellant is 




If Wp is constant then the burning time t^ is 
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Thus any time the condition R - r = O is not satisfied t^ will be changed. 
For the case where r< R the fuel will be exhausted first and 

x f w t fLLL) 

The result is that there will be residual oxidizer on board at burn out. 

(-o = w„-fx p dt = W(ft,-r) 

)c R+l 

For the case where r R, the oxidizer will be exhausted first and there will 
be an excess of fuel on board at burn out. 
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Figure Il-la shows the variation in L as a function of r/R for the assumed 
design given in Appendix I. 

Considering an idealized vertical trajectory where drag has been neglected 
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\jt) » 9 Is f In ~~~T ~ 9^ 
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Substituting t^ for t in the above equations the velocity and altitude at burn 

out can be found as functions of r/R. Figures Q-l b and c shows the performance 

variations as a function of r/R for the assumed design given in Appendix I. 

If some performance level is an absolute requirement it is then necessary 
to carry enough propellants to insure that neither fuel nor oxidizer will have 
been exhausted prior to reaching this desired performance level. This consider- 
ation may add much to the overall size of the rocket due to the high growth 
factor in rocket propulsion systems. 

This is best illustrated by considering the velocity equation in the following 
form 



Where I and are assumed to be constant Exclusive of added tankage 
weight this factor is 13 for the design shown in Appendix I . As a general 
conclusion it is obvious from the above equation that for a lift-off acceleration 
of l/2g this growth factor is in the order of three times the ratio of initial 
weight to burn out weight. 
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For a given burn out velocity requirement it can be shown that 
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FIGURE II- 1 Effect of improper propellant utilization on rocket 
performance at burn out 
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THRUST CHAMBER CONSIDERATIONS 

Thermochemical equations may be used to obtain analytical expressions 
for the conditions existing within the thrust chamber relative to the performance 
of a rocket motor. Given some proposed combination of fuel and oxidizer and 
an assumed chamber pressure it is possible to obtain equilibrium values of 
T.ft, , and l' as a function of mixture ratio. The actual values of the para- 
meters will differ for different propellant combinations but in general are 
observed to follow the general pattern shown in Figure HI-1. Specific examples 
of performance charts for several different propeUant combinations are shown 
in Reference 16. 

The results of the thermochemical analysis will not vary appreciably 
with the chamber pressure assumed providing the combustion equilibrium 
equations are not pressure sensitive. With the assumption of pressure 
insensitivity the chamber conditions can be considered as a function of 
mixture ratio alone. 

One of the basic performance parameters of a rocket motor is the 
characteristic velocity c*, which can be expressed as a function of the com- 
bustion product gas properties as 




With T,fo , and if reduced to functions of mixture only it is seen that c* 
is also a function of mixture ratio. For many bipropellant systems Y may 
be considered constant over the range of r under consideration, thus 
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A plot of fi~ ag a function of r typically results in a double-valued 



cure such as shown in Figure m-2 . From the plot 



o£ '(iF 



derived from 



thermochemical considerations it is possible to obtain an analytic expression 
for c* in terms of mixture ratio by the use of various curve fitting techniques. 



Figure III— 2 shows the variations of 







% 



with r for a liquid oxygen-ethyl alcohol 



(with 25% water) propellant combination. Here a parabolic approximation has 
been assumed to be valid over a 10% range of mixture ratiq as measured from 
X - 1. 37. The resulting values of c* over this band of mixture ratios are 
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shown in Figure H-3. Based on the parabolic approximation for 
given as 

(£-<)- -M(i-i)' 

The dynamics of the combustion chamber are developed from the consider- 
ation of the requirements for mass flow conservation through the chamber. 




Ill -3 



Referring to the above sketch; under steady state conditions the flow through 
the chamber may be expressed as 

(t-t c ^ = rn b (i) - = m 

where t„ represent the time delay between injection and combustion. For non- 
steady conditions the flow is expressed as 

If the time lag t c is neglected and the flow at each point expressed as a 
perturbation about the steady state value the equation becomes 

- 

which can be rewritten as 

_ _d_ 

rr\ dt 

Introducing a new time variable TT such that T = t/t s , then 
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Here t 0 is chosen as the stay time of the flow in the chamber. Therefore 

AlTc — d. ULc + e. 

\r\ ^ c 

The fluid mass in the chamber may be represented as 



m c = +■ &m c 
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where 
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dr 



Thus 
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the chamber is 


given by 



and by substituting the previously determined expression for c* 

- K<r Vc Pc — 



lt\ = XL? -4- = 
c R 1 cT ^ ci 11 



and 



= 2 AC* 



For steady nozzle flow 

ft 
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and 
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AC 



AA . 

m Pc C* 

The flow through the chamber is now expressed in linear form for small 
perturbations about the steady state condition as 

AA c = Ol / hjjc _ Z &C* \ ^ A_fi: _ 
m dr ^ P c c # ) Pc c> 

Using the operator s to signify d/d 7! the flow equation becomes 
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The injector flow may be considered in terms of the individual propellant 
flows and the mixture ratio from the following relationships. 



r = 





hr\ +• tf\r 
o 



Thus the injector flow may be expressed in two forms 



= m* ( i* <r) 



Therefore 
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or in terms of volumetric flow where 

Am i _ 4 

™ Or ) 9<.« 

Substituting this express in the chamber flow equation gives 

AA.* /. A?c 



rt\ = pP 
A Qa 
(1+0 9a 



lWl _ + .-^1 = (i*s)^ - 0+«)iS 
M<?u ^ c 



If the variation in c* is assumed to be uniform throughout the chamber and, 
as before, the combustion lag time is neglected the last term on the right 
side of the above equation may be written in terms of injection mixture ratio 

0 + 2S V )~ = 4 r — (j + *s) Ar 

c* O' ctr ^ 



and in final form the chamber equation becomes 
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u + & r t = K„ &<)*. K„ ^ 



where 
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MIXTURE Ratio / r 



FIGURE ffl-1 Typical performance parameters as a function 
of mixture ratio for a liquid bipropellant rocket 
engine 
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FIGURE III -2 Variation of with mixture ratio 
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FIGURE ID-3 



Variation of characteristic velocity, 
c* with mixture ratio 
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